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Abstract 

Large deviation estimates for the following linear parabolic equation are studied: 
{x)D'^uj + b{x) ■ Du+ ^{u{x + y) - u{x) - {Du{x) ■ y)t{\y\^^ 

iv)] 



where is a Levy measure (which may be singular at the origin). Assuming only that some 
negative exponential integrates with respect to the tail of /x, it is shown that given an initial 
data, solutions defined in a bounded domain converge exponentially fast to the solution of 
the problem defined in the whole space. The exact rate, which depends strongly on the 
decay of ^ at infinity, is also estimated. 



Resume 

Nous donnons ici dcs estimations de grande deviations pour I'equation parabolique lineaire 
suivante: 



(a(x)D'^u^ + b{x) ■Du+ I {(u(a; + y) - u{x) - {Du{x) ■ y)l{\y\<^{y)} d^l{y) , 



du „ 

oil fi est une mesure de Levy (possiblement singuliere a I'origine). En supposant seulement 
qu'une exponentielle negative est integrable par rapport a la queue de distribution de la 
mesure /Lt, nous montrons que, etant fixee une donnee initiale, les solutions definies dans un 
domaine borne convergent vers la solution du probleme dans I'espace tout entier a un taux 
exponentiel. Le taux de convergence exact, qui depend fortement du comportement de /x a 
I'infini, est egalement estime. 
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1. Introduction 



The equations we eonsidcr in this paper take the following form: 

^ = Tr(^aix)D^u^+bix)-Du+J^^ [{u{x+y)-uix)-iDu{x)-y)l{\y\^iy{y)} dfiiy) . (1.1) 

Recently, equations like (jl.ip involving Levy-type non-local terms have been under thorough 
investigation by many authors and in many directions. Indeed, they present challenging 
problems which are not covered by the existing "local" vision of pde's. References are too 
many to cite, those which are more closely related to our concern are given below. 

Our main goal here is to obtain some estimates on how the solutions Ufi that are defined in 
the ball Bj^ = {\x\ < R} approach the solution in as i? — *■ oo, in the spirit of 0]. To 
sum up briefly, the main results of the present paper show that provided /i has at most an 
exponential tail, a bound of the following kind always holds as i? — > oo: 

\u-UR\{x,t) ^ c-«/("^«)(i+°(i». (1.2) 

The behaviour of / depends on the tail of fi, but in any case / cannot grow faster than 
hnearly. This has to be compared to the classical exp(— i?^)-type estimate associated to the 
heat equation 0|: the presence of non-local terms implies an order exp(— i?lni?) at best. 



Typical examples - As for (|l.ip . a first example we have in mind is the convolution 
equation, 

— = J*u-u. (1.3) 

Here a,b — and /x is a Levy measure with probability density J, which is symmetric and 
integrable - see 13, 14 1 . In Q , the authors gave partial answers for this equation when J 
decays strictly faster than an exponential at infinity. Our results here include and generalize 
those given in 0. Be carefuU that writing the equation under the form (|1.3p requires J 
to be of unit mass, and moreover that in this case /x(z) = J(— z), a detail which has to be 
taken into account when dealing with non-symmetric kernels. 

A more sophisticated example is the following: 

du f ( 1 c ^ ^ ^ 

— = Au + b- Du + VN. / \u{x + y,t)~u{x,t)\^^j^dy, 

Ot Jrn I. J \y\ 

where P.V. stands for principal value. This is the same as putting the compensator term 
Du(x)-y since the measure is symmetric. The measure J(y) = e ~'^'/l2/l^^"' called tempered 
a-stable law, appears for instance in finance (20| . 

It is important to notice that our framework (which will be fully established in Section (2]) 
does not allow to treat fractional Laplace-type nonlocal terms for which J{y) = l/|y|^+", 
a e (0, 2) since we require the tail of J to go to zero at least exponentially fast. 



Existence and uniqueness of solutions - Concerning the existence and uniqueness of 
solutions of p.ip . first notice that in the case of constant coefficients a and b, this can 
be derived from a Fourier analysis of the equation. In the absence of differential terms, 
i.e., a = b = and if fi is not singular at the origin, various results about existence and 
uniqueness of solutions in or bounded domains can be found in [^, [lol . 
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Finally, in the case when a(x) = acr* with a and b Lipschitz, existence and uniqueness of 
bounded solutions are proven in for the elliptic (stationnary) version of the equation under 
some assumptions on fi. With little modifications (which essentially amount to changing 
the (— 7w)-term in the equation for a wt-term), the same result holds true for the parabolic 
version, and gives existence and uniqueness of solutions in or in a bounded domain. See 
also P.flll.[l3| for some other results. 

We shall not derive here a full theory of existence, uniqueness and comparison for (jl.ip 
which is not the central point of this paper. Rather, we will simply assume that given an 
initial data uq continuous, positive and bounded, there exists a unique viscosity solution 
u{x,t) of (jl.ip with initial data uq- We also consider for any R> the solution uji of (jl.ip 
in Bji with initial data uq and u{x,t) = for any x > R (see @ for details) and assume 
that for any i? > 0, there exists a unique viscosity solution it^. Starting from this, we will 
now estimate the difference (u — Uf(). 

Large deviations - Estimate (|1.2p may be seen as a large deviations result if one considers 
the probabilistic viewpoint associated to the equation. We shall not enter here into details 
about Levy processes and refer for instance to [1, [l^ for a probability-oriented approach, 
but let us just mention a few facts. 

A Levy process is a stochastic process (Yt)t^o which has stationary and independent incre- 
ments: for any < s < t, {Yt —Yg) only depends on {t — s) and is independent of (Yf — Yg') 
for non-overlapping time intervals. An important fact is that (Yt) is not required to be 
continuous (in time), so that this type of process can model either continuous diffusions or 
jump diffusions, and a mix of both. Thus, typical Levy processes are the "usual" brownian 
motion, the compound Poisson process and the a-stable jump process. 

Now, if u{-, t) represents the density of the law of Yf , then it can be shown, using the Lcvy- 
Khintchine formula, that the characteristic function of Yj, which is the Fourier transform 
of u, takes the following form 

u(^,i)= e*'^(«), 

where the characteristic exponent ip is the sum of of (i) a brownian motion, {ii) a drift and 
(iii) a pure jump process. More precisely, in Fourier variables: 

^(C) = Ae-e + i?-e + / (e'(«-^)-i-i(e-y)i{|,|<ij) d^y), 

where /i is a Levy measure. By taking the inverse Fourier transform of ut = (p(^)u, we 
recover exactly equation (jl.ip : Brownian motions are associated to Laplacian-type diffusion 
terms Tr {a{x)D'^u) while a drift corresponds to the h{x) ■ Z?M-term. In the case of the 
compound Poisson process, 

^(e)= 1 (e'(«^')-l)j(y)dy, 

where J £ L^ is the jump distribution of the process so that the non-local term can be 
re-written as Jh<u, leading to (|1.3p . Finally, in the case of the a-stable law /i(z) = 
we recover the well-known fractional Laplace operator. 

In some sense, estimate (|1.2p measures how many processes have escaped the ball Bji be- 
tween time i = and t ~T (we refer to Q for a precise proof of this assertion in the case 
of the Laplacian), and obtaining an exponential- type estimate is precisely what is called a 
large deviation result - see the introduction of fol, and for more about this. 
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Notations - Throughout the paper we wiU denote by (cj) the canonical basis of R^, while 
V will refer to any unit vector. The euclidean scalar product is denoted as a ■ and we use [ • | 
for the associated norm in R^. We will use ang(a, 6) to denote the angle between the vectors 
a and h (see Definition 15. ip and a — > O"^ is the transposition. The notation B{x,r) stands, 
as is usual, for the ball of radius r > 0, centered at a; G and we also use the simpler 
notation Br = B{0,r). We also denote Du G the gradient of u and D^u G A^Ar(R) the 
Hessian matrix of u. In all the paper, c(J) denotes any constant which only depends on 
J, possibly varying from one line to another. The notation Ji =^ J2 essentially means that 
Ji ^ J2 in R^, but we add some strict separation inside a small region (see Definition l2.4p . 
We denote by a A 6 the infimum of a and b. 

Main results - We have first a general result. Theorem 13. H which gives a theoretical 
bound in terms of a rate function loo, typical of large deviation estimates: as i? ^ 00, 

\u-UR\{x,t) ^ e-«^~(^/^'*/«)+°W«. 

The rate function I^c satisfies the limit Hamilton-Jacobi problem 

dJoo + H{DIoo) =0 inBiX (0, 00), 

/oo=0 on 9Bi X (0,00), (1.4) 

/oo(2:, 0) = +00 in _Bi, 

where the Hamiltonian associated to (|1.1[) is given by: 

Hip):=Tr{Ap®p) + B-p+ f (cPy ~ 1 ~ {p ■ y)ti\y\^,}iy)) dfi{y) . (1.5) 

The idea of the proof is to suitably rescale the problem for u—ur (which is not the same as in 
the case of the local heat equation, though) and to derive an equation for the log-transform; 
then pass to the limit as i? — > cxd. 

The Lagrangian L associated to H is defined as usual as 

L{q) := sup {p-q-H{p)], (1.6) 

peR" 

and using a Lax-Oleinik formula, see [l9| . we obtain a semi-explicit expression for loo, 
see p.2p in terms of L. 

In our present framework, obtaining (jl.4p is more involved than in since we have to deal 
with the differential terms and the singular part in H . Hence, we face several difficulties 
that come from the generalization of the equation, which imply non trivial adaptations and 
new developments of some aspects of the theory. Let us briefly explain the main points we 
have to deal with: 

(i) Singular measures at the origin: using kernels with singularities requires first a suitable 
concept of solution; we use here the notion of viscosity solutions derived in for 
Levy-type operators, which allows us to handle this situation. Then we introduce 
what we call the essential Hamiltonian, H^^^ , wiping out the singularity at the origin, 
see Definition 12.31 and prove that asymptotically H and H'^^^ are equivalent (and the 
same holds for the associated Lagrangians) . 



4 



(ii) Non symmetric kernels: wc want to estimate L in order to obtain an expression for 
lao ■ This needs an extra effort in this case: since we do not have symmetry, the vectors 
p and DH{p) do not necessarily point in the same direction. Hence dot products can 
not be simphfied and wc have to carefully analyze the angle between vectors. Wc prove 
some results valid for non symmetric kernels J at the logarithm level, and others using 
the smallest symmetric kernel above J. 
{in) Differential terms: the first two terms in equation p.ip do not introduce much diffi- 
culties with respect to Q but we have to deal with them with the notion of viscosity 
solution. To be able to pass to the limit, we will assume that a and b have limits at 
infinity. 

(iv) Possibly infinite Hamiltonians: in the case when J{y) ~ c~"'^' as \y\ — > oo, H 
is infinite outside B^, which makes the analysis of \1A\ much more delicate, since 
an initial layer appears. Thus, one of the major contributions of this paper is to 
provide existence, uniqueness and comparison results for the Hamilton-Jacobi equation 
ut + H{Du) = when H its infinite outside the a ball. 

Organization - After a preliminary section. Section^ in which some properties of the gen- 
eralized non-local equation as well as for the Hamiltonian are stated, we devote Section [3] to 
the theoretical behaviour of the problem, where we prove Theorem l3.1l Then we concentrate 
our efforts in finding a behaviour for L, that allows us to give more explicit convergence 
rates for u — ur. We deal with this issue in sections 21 and [SI where we consider different 
types of measures (compactly supported, intermediate, exponentially decaying), which lead 
to different behaviours. Section [S] especially deals with critical kernels for which H is not 
finite everywhere; so we devote a big part of this section to treat the Hamilton-Jacobi equa- 
tion with such H . Finally in Section[7]we collect some further results: we give some explicit 
calculations of rates and consider totally asymmetric measures in 1-D. We also explain how 
our methods allow to obtain some estimates for the related nonlocal KPP-cquation. 



2. Preliminaries 

Let us focus now on the concrete hypothesis that concern the equations we have in mind. 
For jjL, a, b, A, B in p.ip and p.Sp there are essentially three assumptions we shall make 
throughout the paper: 

Hypothesis 1. The measure is a Levy measure with density J(-) satisfying 

' J : \ {0} ^ R+ is C^-smooth, 

< 3/?o>0, supp(J) D Bpo, ^2.1) 

/ (lA|ynj(y)dy<oo. 

Hence in the sequel we rewrite (|1.5p as 

H{p):=Tt{Ap®p) + B-p+ f (0^^ - 1 - (p • 2/)l{|,|<i}(y)) J(y) dy . (2.2) 

The last assumption in (|2.ip implies in particular that the integral in (|2.2p over {|y| < 1} 
converges. 

Hypothesis 2. The N x iV-matrix a(-) is continuous and nonnegative, &(•) is a continuous 
vector field in and we also assume that the following limits are well-defined: 

A= lim a(.T) e 7W„(R), B= lim b{x) eS.'^ , (2.3) 
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(notice that necessarily A is nonnegativc) . 

Except in Section [SI we also make the following assumption on J: 
Hypothesis 3. 

V/3>0, / e^^y^J{y)dy <oo. (2.4) 

J{\y\>po/2} 

This ensures that the exponential part of the Hamiltonian in ()1.5p (or (12. 2p ) converges for 
any p G R^. In Section [6l we shall only assume that 

sup{/3>0: / e''l^lj(2/)d2/ < oo} =/?o > 0, (2.5) 

J{\y\>po/-2} 

which is the case for exponential- type tails, and implies that the associated Hamiltonian is 
not everywhere-defined. 

Most of our results would be also valid in a more general setting, for instance we could relax 
the condition on J, but this is not a major matter for this paper. There are also some 
specific other assumptions that we shall make in Section [5l 



2.1. Viscosity solutions of the equation 

Equations like (|l.ip can be treated by various ways, depending on the coefficients a(-) and 
6(-) and whether the measure fi is singular at the origin or not. Since we make only general 
assumptions, the best tool in order to deal with differential terms and a singular measure 
is provided by the notion of viscosity solutions. 

Essentially, the notion of viscosity solution allows to give sense to the differential terms 
as well as the integral term without knowing a priori that the solution is regular, using 
comparison with smooth test functions. Notice however that in the integral term, one only 
needs to replace ti by a smooth function in a small neighborhood oi y — 0. We refer to 
BE I 121 for more results on viscosity solutions in presence of singular measures. So, if we 
set 

Eq[w, 0, (S](a;o,to) := ^{xo,to) - Tr (^a{x)D^(f>{xo,ta)j - b{xo) ■ D(j>{xo,to) 

- / liHxa + y,to) - (t)(xo,to) - (D(t){xo,to) ■ y)\j{y)dy 
J{\v\<s} ^ 

- / \{u{xo + y,to)-u{xo,to)-{D(t){xo,to)-y)t{\y\<i}\j{y)dy, 
we then introduce the following definition: 



Definition 2.1. A locally bounded u.s.c. function u : x [0, oo) ^ R is a subsolution of 
iff, for any 5 G (0, 1) and any test function £ C^(M^ x M+), at each maximum point 
(xo, io) £ X 1R+ of u — 0, we have 

Eq[ii, 0, (5](a;o,io) ^ 0. 

A locally bounded l.s.c. viscosity supersolution is defined in the same way with reversed 
inequality and min instead of max. Finally, a viscosity solution is a function which the 
upper and lower semicontinuous envelopes are respectively sub- and supersolution of the 
problem. 
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When intial and/or boundary data are involved, the definition takes into account that at a 
boundary point, cither the inequation has to hold, or the boundary data has to be taken in 
the sub/super solution sense. If CI is the domain of definition, / is the boundary data (in 
the sense of nonlocal equations) and u{x,0) = uo{x) is the initial data of the problem we 
define: 

Definition 2.2. A locally bounded u.s.c. function u : M.^ x [0, oo) — > R is a subsolution 
of with boundary data u = / in (R^ \ $7) x (0, oo) and initial data u{x,0) = uo{x) 

iff, for any (5 > and any test function cj) G C^(R^ x (0, oo)), at each maximum point 
(xo, to) & CI X [0, oo) of u — (/), we have 

(xo^to) en X {0,co) => Eq[u,(/),(5](xo,to) < 
(xo^to) e dCl X {0,oo) min (^Eq[w, 0, 5](xo,to) ; u(xo,<o) - /(a;o,io)) < 

{xo,to) e CI X {0} => u{xo,to) sS,uo{xo) 

A supersolution is defined with reversed inequalities and min/max changed accordingly, and 
a solution is such that its l.s.c. /u.s.c. envelopes are sub/super solutions of the equation. 

2. 2. Hamiltonians 

We shall show now some properties of the specific Hamiltonians we consider: 



Lemma 2.1. Let J be a kernel satisfying ()2.ip . Then the Hamiltonian H is superlinear and 
strictly convex, and thus the associated Lagrangian L is well-defined, convex and superlinear. 



Proof. A straightforward calculus shows that in the sense of matrices, 

D''H{p) [ (y ® y)eP-yj{y) dy > 

SO that H is strictly convex. Now, one can easily check that 

H{p) = 0{\pf)+ f efV(y)dy. 

Indeed, the differential terms are clearly of the order of |pp and \p\ respectively, while the 
integral over |?/| < 1 can be bounded by: 

f (ePy-l-ip-y))j{y)dy^C\p\^[ + 0{p ■ y))J{y) dy ^ C\p\\ 

J{\y\<i} ^ ^ h\y\<^ ^ 

On the other hand, 

eP V(2/)d2/^ \ eP yjiy)dy;? c{J)cP»\P\/^ , 

{|y|>i} J{p-y>po\p\/2} 



so that indeed H is superlinear. It is well-known (see [2l| for instance) that if H is strictly 
convex and superlinear, so does L. □ 

We shall go a bit further in this direction, using that the main contribution of the Hamilto- 
nian comes from the exponential term in the integral. So let us define the essential part of 
the Hamiltonian and the corresponding Lagrangian: 



7 



Definition 2.3. We denote by 

:= /" cP '"J{y)dy and L''''''{q) := sup {p ■ q - H''''''{p)} , 

J{\y\>P(,/2} pSK" 

where po is defined in (|2.ip . 

Tlie reason wliy we integrate over {|y| > po/2} is that we want to avoid the singularity at the 
origin, but we also want to be sure to integrate within the support of J when it is compactly 
supported. The condition supp (J) D Bp^ ensures that {|y| > po/2} n supp (J) ^ 0. 

We have first a basic estimate, similar to the one that was used in Lemma |2. II 

Lemma 2.2. Let J he a kernel satisfying (j2.ip and (|2.5p . Then we have, 

H'^'^'^ip) ^ c(J)e'"'lfl/2, 
p-DH-'^{p) ^ c(J)(po|p|/2)e''°lfl/2 + 0(|p|). 

Moreover there exists e > such that for any unit vector v, 

^T^2^css(^)^ ^ ec( J) e IpI'"'/^ 

and the same result holds for D'^H = A + D'^H'''''' ^ D^iJ^^s 
Proof. The first estimate is easily obtained: 

H'^^^ip) ^ [ e^yj{y)dy 

J {p-y>po\p\/2} 

^ c(J)e''°IPl/2. 

Notice that the constant c( J) may be small, but it is still positive since J is continuous and 
positive inside {p ■ y > pq\p\/2\. 

For the second estimate, the proof is essentially the same: 

p-DH'^%p) ^ f {p-y)eP yj{y)dy+ f {p ■ y) y J{y) dy 

J{p-y>po\p\/2} •'{p-y<o} 

^ ciJ)ipo\p\/2)e'"'\p\/^ + Oi\p\) - c(J)b| , 
where in this case c(J) = I{\y\>i} \y\Jiy)'^y < °° from (|2.5p . 

Finally, for any unit vector of the canonical basis of M^, we compute with the same 
decomposition: 

efD^H''%p)e,= [ e^' V(y) dy [ \y,\^ eP-y J{y) dy . 

J{\y\>po/2} J{p-y>Po\p\/2} 

We notice that the set {p ■ y > pa\p\/2} contains the truncated cone: 

f p-y .3 I , 3po \ 

hence for some e > small enough (which can be chosen independently of e^), 

B:={p-y> po\p\/2}n{\y,\^^e}^ill, 
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and since B fl Bp^ ^ 0, the integral of J on B is positive. So we can estimate the integral 
from below as follows: 

eJD^H''^%p)e^ ^ ee^P^P"^^ [ J{y)dy 

Jb 

^ ec{J)c^P^P°'^ . 

So the result also holds for any unit vector v. □ 
The next estimate will be essential in the sequel: 

Lemma 2.3. Let J he a kernel satisfying (|2.ip and (|2.4p . Then, for any 7 e (O,po)j 

F-»^^ + ^:^^+e^H. (2.6) 

7 l\p\ 

As a consequence, H°^^{p) is negligible in front of p ■ DH°^^{p) as \p\ 00, 

Hm ^i-i— . 

IpHoo p ■ DH''''^{p) 

Proof. We begin with splitting the integral in two terms, 

H-''{p)= [ eP yj{y)dy+ [ qP ^ J [y) Ay = h + h 

The first integral is estimated as follows: 

/i < / e^lPlj(2/)d2/^c(J)e^lfL 

J{p-y<7|p|}n{|y|>po/2} 

Now, notice that since 7 < po, even if J is compactly supported, we are sure that the integral 
I2 concerns a region where J is not zero, so that we can indeed control it by p • DH^^^ij)): 



I2 f {p-y)ePyj{y)dy 

= ^ / {p.y)ePyj{y)dy-^ f 



{p-y)ePyj{y)dy 

{p-y^ilpl} 

^ p-i^g-(p) _ 1 f ^p.y~^^P yj^y)^y_±_ j ■ ^ ) C ^ ^ ^ ( ^ ) d ^ 

7bl 1\P\ J{p y^O} 

^ p ■ DH-^^p) ^ cjJ) 
l\p\ 1 

Summing up Ii and I2 gives the first result of the lemma. 

For the second part we first fix /i > 1, 7 = and choose \p\ > /-t/po so that 7 e (0,po)- 

Using estimate (|2.6p . we obtain: 

i^-»<^ + ^:^^ + e^ 
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which implies 

p-DH'^'^p) c{J)\p\ ^ic^ 



Since H°^^{p) — > oo superhnearly (see Lemma [2?2|) we then have 

. pJ^IT^ 
lim ini — ^ u . 

But since /i > 1 is arbitrary, we conclude that the limit is +00, hence the result. □ 

Remark 2.4. If DH°^^{p) would grow faster than an exponential, estimate (|2.6p would be 
sufficient to conclude that H'^^^{p) is negligible in front of p ■ DH''^^{p). Actually this will 
be the case when J is not compactly supported, see Lemma 15.61 



The following lemma proves that essentially, the estimates we will produce in this paper do 
not depend on A, B and neither on the behaviour of J near the origin. 

Lemma 2.5. For any J satisfying (|2.ip and (j2.4p . as \p\ — > 00, we have 

H{p) - , DH{p) - DH''''^p) and L(p) ~ . 

Moreover, if J is symmetric, then iJ^^*^ and L°^^ are also symmetric. 

Proof. We split the Hamiltonian as follows: 

H{p) = Tr {Ap ® p) + B ■ p + H'^'^^p) + [ ( ^ - 1 - {p ■ y)) J{y) dy 

J{\y\^po/2} ^ ^ 

- / (l + ip-y))j{y)dy- [ J{y)dy 

= H'''%p) + Tt{Ap®p)+B-p+ I (^l:^^i + 0{p-y)))j{y)dy + c{J) 

J{\y\^po/2] ^ ^ ' 

^H^^%p) + 0{\p 
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Since the Hamiltonians behave at least exponentially, Lemma [2.21 the 0(|pp) is negligible 
by far, and H and H'^^^ are equivalent. The calculations for DH and DH'^^^ are similar: we 
get 

\DH{p)-DH''^p)\=0{\p\^), 

while \DH{p)\ grows at least exponentially. Notice that since H and H'^^^ are equivalent. 
Lemma implies that both are negligible in front of DH and DH'^^^. 

More involved is the proof that L^^^{p) ^ L{p). Let us notice first that the point pq = po((j) 
where L(q) reaches its sup goes to infinity as \q\ — > 00. Indeed, this comes from the fact 
that q = DH{po) and that DH{p) grows at least exponentially, see Lemma [2.21 The same 
holds for DH''''\ 

Then we proceed as follows: for \q\ large enough, since the sup below are attained for \p\ 
also large, there is some constant C > such that 

sup {p-y~ H''%p) - C\p\^} ^ L{q) sup {p • y - H^'^p) + C\p\^} . (2.7) 
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Hence if we denote by pi ~ Pi{q) the point where the sup on the right is attained, we have 
q = DH'^''''{pi) + 'SC\pi\pi, while q = DH{po) for the sup in L{q). Now, 

q = DH-'%p,) + 3C\p,\p, = DH-'%p,) + 0(|pi|2) = DH{p,) + 0{\p,\^) = DHipo) , 

and from this we wiU concKide that pi (q) ^ pq (q) . More precisely, 

0(|piP) = DH{p,) - DH{po) = D'HiO ■ [Pi -Po) , 

for some ^ in the segment [po,Pi]- We take the scalar product with [pi — po) and use 
Lemma [2.21 to get an estimate from above: 



bi -por(£c(J)el«l''«/2) < [D^HiO ■ (pi -Po) ■ (Pi-Po)) ^ C\pi\'\pi - po\ 

From this follows that for some constant still noted C > 0, 

e-l«IPo/2 

\Pl -Po\ ^ C ■ f2 . 

Noticing that |^| ^ minjlpol, \piW ~^ oo, this estimate implies that as j^l — > oo, 

(J g -inin{|poMpi|}po/2 



min{|po|, bill bip min{bo|, bill 
Hence we have proved that pi{q) ~ Poiq)- 

The same calculation is valid for the sup on the left of (|2.7p . so that we conclude that 

L-'^iq) = pi • DH^'^pi) - H-'%pi) - PO ■ DH{p„) - H{po) = L{q) . 
About the symmetry see [^, Lemma 2.4]. □ 

We finally prove a result which allows us to compare the Lagrangians when the kernels 
are ordered. But for some technical reason, we need the kernels to be strictly ordered in 
Bp \ -Bpo/2, hence we introduce the following notation: 

Definition 2.4. We say that Ji and J2 are essentially ordered and we denote Ji =^ J2 if 
there exist a, 6 > such that po/2 < a < b < po and 

Ji < J2 in , My) < My) in {a < \y\ < b\} . 

Lemma 2.6. Let us assume that Ji,J2 satisfy (|2.ip and (j2.4p and that Ji =^ J2. Denote 
by Hi, H2, Li, L2 the associated Hamiltonians and Lagrangians. Then for p,q big enough: 

Hi{p) H2{p) and Li{q) ^ ^2(9) • 

Proof. The inequality Hi^^{p) ^ H2^^{p) is always true, since the integrand in H''^^ is 
always positive and Ji ^ J2. But as such, this is not enough to derive the same inequality 
between Hi and H2. To this end, let us notice that since Ji A J2, then there exists e > 
and pq/2 < a < b < Po such that 

My) ^ My) - exiy) in M^, 

where x(^) = ((^ ~ l-'^Dd^^l ^ ■ 
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Then we shall estimate the difference H^^^ — i/l^** as in Lemma [121 
Hrip)-Hrip) ^ -s [ cf «x(y)dy 

J{a<\y\<b} 

^ -£ / eP-yx{y)dy 

J {a<\y\<b}n{p-y>a\p\} 

where c(x) is the integral of x over {a < \y\ < b} Ci {p ■ y > a\p\}. It is clear that this set 
has non-empty interior and since x > in this set, then c(x) is a stictly positive constant. 
This proves that not only are Hf^^ and H2^^ ordered, but moreover the difference between 
the two is at least of the order of an exponential. 

Then we use, as in the proof of Lemma [2.51 that Hi{p) = Hf^^{p) + 0{\p\^) and the same 
for H2, so that 

H,{p) - H2{p) = nrip) - Hrip) + 0{\p\^) ^ -£c(x) e I + 0{\p\^) . (2.8) 
Thus for \p\ big enough, wc have indeed Hi{p) ^ H2{p). 
The inequality for Li and L2 follows by definition. □ 

Remark 2.7. [i) This last result is also valid if Ji satisfies (|2.5p with J2 satisfies (j2.5p 
with (32 ^ /?! and wc consider \p\ < (32- (ii) Notice that if wc only assume Ji = J2 outside 
Bp and Ji < J2 in Bp, then Ji =<! J2. In this case, the tails of both kernels and their support 
(if compactly supported) remain exactly the same so that essentially, all the estimates we 
give in the sequel are the same for both kernels. 

2.3. Viscosity solutions of the limit Hamilton-, J acohi equation 

In Section [31 wc will have to study the following Hamilton- Jacobi equation with Cauchy- 
Dirichlct boundary values, 

dtl^ + H{DI^) ^0 inBix(0,oo), 

7-^ = on OBi X (0,cx)), (2.9) 

where we assume here that H e C(M^;]R). We refer to Section [6l for the case when the 
Hamiltonian could take infinite values. Let us first recall the definition of viscosity solutions 
for this equation (sec for instance [13] )■ 

Definition 2.5. A locally bounded u.s.c function u : Bi x [0, 00) — > R is a viscosity subso- 
lution of (|2.9p if for any C^-smooth function (p, and any point (xo,to) & Bi x [0,oo) where 
u ~ ip reaches a maximum, there holds, 

xqEBi =^ dt(p + H {Dip) s^O at {xo, to), 
xq e dBi => min {dt(p + H{Dip) ; w} ^ at (a;o, to) , 
io = => u ^ A in i?i . 

A locally bounded l.s.c. function is a viscosity supersolution if the same holds with reversed 
inequalities and min replaced by max at the boundary. Finally a viscosity solution is a 
locally bounded function u such that its u.s.c. and l.s.c. envelopes are respectively sub- and 
super-solutions of (|2.9p . 
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Notice that in general, the initial data as to be understood in the relaxed viscosity way 
(either the data is taken, or the sub/super solution condition holds) but it is well-known 
(see for instance 0] ) that in the case of continuous Hamiltonian, this condition is equivalent 
to the one we use here. Keep in mind, however, that this will not be the case in Section [6l 

Since H is convex we have the following representation: 

Lemma 2.8. // J satisfies (|2.1[) and (j2.4[) . then there exists a unique viscosity solution 
of (|2.9p . which is given by 

I^{x,t) = A A mm \tL(^^—^)\ in Bi x (0, oo) . 

Proof. The assumptions made on J imply that both H and L are finite everywhere, convex 
and super- linear so that uniqueness holds for this problem, see [Toj . Then the Lax-Oleinik 
formula in the bounded domain Bi x (0, oo) gives: 

/^(a;,i) = min (tif^— + a| A min ( {t - s)l(^^—^)] . 

^ ' |y|^l I \ t J i (!/,s)eai3ix(0,t) I Vt-s/i 

Using that L ^ and L{Q) = 0, we obtain that the first min equals A. Then using that the 
function r L{cr)/r is increasing, the second minimum is attained for s = so that the 
result holds. □ 

3. Theoretical Behaviour 

The main goal of this section is to derive a theoretical bound, in terms of the Lagrangian 
L, for the error made when approximating the solution, u, of (jl.ip by solutions, wj^, of the 
Dirichlet problem. This result will allow us to derive explicit rates of convergence provided 
we know the behaviour at infinity of L{q) (this will be the aim of the next sections): 

Theorem 3.1. If J satisfies (|2.ip and (j2.4[) . then for any fixed X G and t > 0, as 

R ^ oo, there holds 

\u~UR\{x,t) ^ c-Ri^i.-/R,t/R)+oii)R^ (3^^) 
where the rate function is given by 

I^{x,t) = mm tL(^^-—^) . (3.2) 

yedBi \ t / 

Moreover, for any T > the o(l) is uniform in sets of the form 

{\x/R\ 1 , s$ t/i? ^ T}. 

We shall dedicate the rest of this section to the proof of this theorem. Notice that if one 
only assumes (|2.5p instead of (|2.4p . the theorem still remains valid, but the proof requires 
some finer arguments since the Hamiltonian can be infinite in some regions of . In such 
a case, an initial layer appears that we analyse in detail in Section [6l 
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3.1. Formal convergence 

Let us denote by vr the solution of (jl.ip in Br x (0, oo) with initial value vr{x^ 0) = for 
X G Br and "boundary data" vr = ||mo|1oo for 5^ R. By comparison, ^ u — ur ^ vr 
so that we need only to estimate vr. 

We first rescale the equation both in x and t as follows: 

wr{x, t) = vr{Rx, Rt) for X- e Bi, t ^ . 

Then wr satisfies a rescaled equation in the fixed ball Bi in the sense of viscosity: 

dtWR{x, t) = -1 Tr (a{Rx)D'^WR{x, i)) + b{Rx) ■ Dwr{x, t) 

+ R \wr{x + y/R,t) - WR{x,t) - {Dwr{x) ■ y/R)l[\yi^iy{y)\j{y)dy. 

In order to estimate wr, we follow 0] and perform the "usual" logarithmic transform, but 
we have to rescale accordingly, dividing by R (and not R^ as it is the case for the heat 
equation). So, remembering that for t > 0, wr > 0, let us define 

lR{x,t) = -^ln{wR{x,t)) . 

Then 

dtWR{x,t) = -i?e-^^«(^'*)9t/fl(x,t) , DwRXx,t) = -Re-^^"^""''^ DlR{x,t) , 
D^WRix,t) = -i?e-^^«("'*)D2j^(^^^) ^ R^ e-^^''^'=-^\DlR(g> DlR,){x,t) , 

and 

^^WR{x + y/R,t) -WR{x,t) - {Dwr{x) ■ y / R)t {\y\<i}{y)^ J {y) dy 

^RI^i.+y/R,t) _ ^^RiM.,t) + (i?/^(x) .y)e-^^-(-'*)l{|,|<i}(y)}j(y)dy 

= /■ |e«{-^«(-+^/«.*)+^«(-'*)}-l + (i?/^(x).2/)l^|,|<i}(y)|j(j/)dy. 

As for the differential terms, we have 

- e {^^'^ (a{Rx)D^WR{x, t)^ + ^b{Rx) ■ Dwr{x, i)) 

= Tr {a{Rx)D^lR{x, t)^ - Tr {a{Rx)DlR{x, t) ® DIr{x, t)) - b{Rx) ■ DIr{x, t) . 

We thus arrive at the following equation for Ir (to be interpreted in the sense of viscosity, 
with test functions): 

dtlR{x,t)+ f |e-^{^«(^+^/^^*)-^-(^'*'} -l + (Z?/^(x).y)l^l,l<i}(y)|j(j/)dy 

= - Tr {a{Rx)D^lR{x, t)) - Tr {a{Rx)DlR{x, t) ® DIr{x, t)j - b(Rx) ■ DIr{x, t) 

Using the limits defined in (j2.3|) . the equation formally converges to the Hamilton- Jacobi 
equation dtl + H{DI) = with 

H{p) = Ty[Ap®p)+B-p + y (cf^-l-(p.2/)l||j,|<i})j(y)dy. (3.4) 

In the following subsection we then justify the convergence of Ir towards the solution in 
the sense of viscosity solutions. 
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3.2. Passing to the limit in the viscous sense 

A first problem comes from tlie fact that if wji approaches zero, then may not remain 
bounded. Hence to avoid upper estimates for 7^, we use the same trick as in which 
consists in modifying Iji a httle bit. For any A > 0, let 

Ir{x, t) = -Jl Hwr{x, t) + e -«'^), 

which is bounded from above by A. Let us notice that since equation (jl.ip is invariant 
under addition of constants, /j| satisfies the same equation as Ir, that is, equation (|3.3|) . 



Proposition 3.2. The sequence (Ip ) converges locally uniformly in Bi x [0, oo) as R +oo 
towards the unique viscosity solution of p.9p . 



Proof. We introduce the half-relaxed limits, for x E Bi, t ^ 0: 

I^{x,t) ■.= limsup*I^{x,t) ^ limsup I^{x',t') 

R—tOD (x' ,t') — (x,t) 

and 

/^(x,t) liminf */^(a:,t) = liminf /j|(a;',t'), 

and we shall prove that they are respectively viscosity sub- and super-solutions of the limit 
problem (|2.9p . Then a uniqueness result will allow us to conclude. 

^ 

Let us take d £ (0, 1) and a test function such that / — has a maximum at (xq, to). Up 
to a standard modification of 0, we can assume the maximum is strict so that there exist 
sequences i?„ — + -|~oo and (a;„,t„) (.To,to) such that 

I^^ — (j) has a strict maximum at (a;„, t„) . 

Case 1: the point (a;o,io) is inside Bi x (0, oo). Then for n big enough, all the points 
{xmtn) are also inside Bi x (0, oo) so we may use the equation for at those points and 
pass to the limit. 

We first write down the viscosity inequality for 

^{Xn,tn) ^ —Tv (a{RnXn)D'^ (j){Xn,tn) 

at Hn ^ 

~TT(a{R„X„)D(l){x„,tn)®D(j){Xn,tn)) - b{R„X„) ■ D(j}{Xn,t„) 



{\v\<&} 



I g-i?{0(..+y/fl,U)-0(x„.t„)} _ 1 _ (i?0(z„,t„) • y)} J(y) dy 

|e-«ui(-+^/«.*")-^l(-"^*")} - 1 - {Dct.{xr.M-y)Hy\<i}]jiy)'^y 

'{\y\>i} ^ 
= Difl[(n) +Inti(n, (5) + Int2(n,(5) 

where DifF(n) represents the differential terms, Inti(?T,, 5) is the integral over {\y\ < S} and 
Int2(ri, S) is the integral over {\y\ ^ 6}. 



Let us first remark that passing to the limit in the differential terms is easy, using p. 31) : 
Difr(7i) ^ - Tr (a • D(l){xo, to) ® L'0(xo, to)) - B ■ D(j){xo,ta) . 
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For n big enough, the first integral term can be controled by: 

|Inti(5,n)K pVllL-(Bi(=^o,to)) / |2;|V(2/)d2/^0 asJ^O. 

-'{l-yK'5} 

It remains to pass to the Umit in Int2. To this end we use the fact that, since we have a 
maximum point, for any z G R^, 

Then, we fix e > 0, choose some M > 1 and spht Int2 into two terms as follows: 
Int2(<5, n)^~f |e _ i _ p0(a;„, i„) . y)l{|,|<ij} J(y) dy 

\wR,^{x,i +y/Rn,t„) - WR,^{xn,tn)\j{y)dy . 
\y\»M ^ J 

Since wr is bounded by l|uol|oo, we can choose M big enough so that the second term is 
less than e, independently of n. 

Then we write a Taylor expansion for near point x„: there exists a ^„ S Bm such that 
Int2 - / I e -DH^^,t„yy+TtiD'-HU)vv) _ i _ (l?0(a;„, i„) . y)l^^^^^Aj{y) dy + e . 

J5(i\y\<M^ ' 

Since ^„ remains in Bm and (j) is smooth we have that D^(j)(S^n) remains bounded. Hence, 
we can pass to the limit as n ^ +oo: 

limsuplnt2((5, n) ^ - f f c -DH-o,toyy _ i _ {^D<t>{x^M) ■ v)t {\y\<i]] J {v) dy + £ . 

n— oo J5i^\y\<M ^ ' 

Summing up the various terms, we obtain that for any 5 > Q and any e > 0, there exists 
M = M (e) > 1 such that 

^(xcto) <-Tr(A-I?0(xo,to)®i?0(xo,to)) ~ B ■ D(j){xo,to) 

- [ |e-^*(--*°)-^-l-(i?(/.(xo,<o)-y)l{|y|<i}}j(y)dj/ + £ + 05(l) 

where os(l) represent a quantity that goes to zero as (5 — > 0. 

It only remains to pass to the limit as e,S 0. Since J satisfies (|2.ip and (|2.4p . the integral 
over converges, and we can send M to +oo and obtain in the limit: 

^(a;o, to) ^ - Tr • D^{xo,ta) ® Z?(/.(xo, to)) - B ■ D^{xo, to) 

_ f |e-^^(-«'*")2'-l-(i?0(xo,to)-y)l{|,|<i}}j(y)dy, 



which shows that / at (I'o, to) is a subsolution in the sense of viscosity. 

Case 2: the point (xo, to) is located at the boundary, xq g dBi,to > 0. Then the sequence 
{xn, tn) may have points a;„ either inside Bi x (0, oo), or at the boundary, or even outside 
Bi. If Xn (z Bi, we use the equation as in the previous case while if a;„ G dBi, we use the 
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relaxed boundary condition in Definition 12.21 Finally, if |a;„| > 1, then Iii^(xn,t„) = so 
that in any case, one has 

min {dt(t> + H{D4>) ; I^^ } ^ at (a;„, f„) . 



We then pass to the limit as n ^ +oo and get the relaxed condition for / at the boundary. 

Case 3: the point (a;o,to) is located at to = 0, € The same as in case 2 happens: 
we have either t„ = and then we use the initial condition, or f„ > in which case we use 
the equation. In any case we get 

min {dt4) + H{D4,) ; /^^ - ^} at (a;„, t„) , 

which gives in the limit 

m:uY{dtct) + H{D(j)); 7^ - A} at (xo,to)- 

Now, it is well-known (sec for instance 0, Thm 4.7]) that in this case, the initial condition 


is equivalent to limt^o I ^ A. Actually this can be proved as in Proposition 16.31 using 
that in the present situation, the Hamiltonian is finite everywhere. 

Conclusion: First, the supersolution conditions for are obtained by the same method, 
with reversed inequalities. Then using comparison between u.s.c./l.s.c. sub/super solutions 
for (|2.9p . we get the inequality ^ I , which implies equality of both functions. Hence, all 
the sequence converges uniformly in Bi x [0, T] for all T > to the unique solution I^. □ 

3.3. Proof of Theorem\^ 

This result only comes from the fact that for any A > 0, by construction 

7^^ = inf(7, A), = inf(/. A), 

with 

I{x,t) := limsup t), Lix,t) := liminf */i^(a;, i). 

The fact that ~ together with Lemma [2.81 vields the result for fixed {x,t), passing 
to the limit as ^ — > cx). 

Now, the convergence of Ir to loo is locally uniform in Bi x [0, oo) so that for any T > 
as long as x/i? < 1 and ^ t/R ^ T, 

\lR{x/R,t/R) - Ioo{x/R,t/R)\^ o{l) ^ as i? -> cx) , 

where o(l) is uniform with respect to x and t as above. Thus estimate indeed holds, 
which ends the proof. Notice that at t = 0, both Ir and loc arc infinite (which corresponds 
to A = +oo), but anyway, the difference Ir — loo remains uniformly controlled. 

Remark 3.3. At this stage we have an estimate valid up to the boundary of Br. In the 
next sections we shall derive more explicit estimates only for |a;| ^ OR, with 9 E (0,1), 
because we use the asymptotic behaviour ot L{q) as |g| —^oo. 
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4. Compactly supported kernels 



In this section we prove that a general "i?lni?" bound is valid for compactly supported 
kernels, extending the symmetric and regular case proved in In order to take into 
account the possible asymmetry of the kernel, we define below for any unit vector the 
size of the support of J in the direction v. 



Definition 4.1. For any e with \u\ = 1, let 

p{i^) := sup {r > : J (ru) > 0} . (4.1) 



Notice that since J is continuous, for any r > close enough to p{iy), J is positive in a 
neighborhood of rv. Notice also that if J is symmetric, then p^iy) = p, the radius of the 
support of J. We shall first derive a bound from below for non-symmetric kernels in the 
logarithmic scale: 



Lemma 4.1. Let J be a continuous compactly supported kernel, let v = and define 

p(y) as above. Then we have: 

ln(w- Dff , , 
lim inf , , ^ 1 . 4.2) 

\p\^oc p{iy)\p\ 



Proof. Let us first choose a unit vector i', define p{i') by (|4.1[) and consider p G M.^ going 
to infinity in this direction: p/\p\ = v and \p\ — * oo. We begin with writing 

p-DH''%p)= f p-yeP yj{y)dy+ [ p ■ y y J {y) dy 

J{p-y^O}n{|y|>po/2} J {p-y>0}n{\y\>po/2} 

The first term is bounded by from below by ~c\p\ and for the second one we define for e > 
and /3 < 1, the set 



CIp = {2/ : ^ ^ Ii/K p{^). p-y^a- e)\p\\y\ ^ 0} n {\y\ > po/2} . 

Hence 

/ p.yePVJ(y)dy^\p\^{l-e)e^'-'^\P\'^ f J{y)dy 

J{p-y>0}n{\y\>pa/2} P J C^i, 

>C(e,/3)W^(l-e)e(-^)H^. 



Notice that C{e,/3) > since J is continuous and positive near p{v)i'^ even if this constant 
could be small. Summing up, 

p.DH'^^\p)^ C{e,l3)\p\^{\-e)e'^^-'^\p\'i^ -\p\c^ KC{e,mp\^{^^ 
for some constant K. Therefore, we obtain for every (5 and e 

^. . \n{p.DH°^%p)) . /lnC(e,/3) ln(^^|p|(l - s)) Eif\p\{l - e)\ i-e 

lim mf , , , , > hm mf ' i i i — 



IpHoo p{v)\p\ IpHoo \ p{v)\p\ P{v)\p\ P{i^)\p\ / P 
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Now, letting e and /3 ^ 1 wc conclude that (|4.2p holds. □ 
In other words, using Lemma [^751 and [^751 we have obtained 

lim inf 1 . 

kHoo p[u)\po{q)\ 

Then in order to have a more explicit bound using Theorem 13.11 we shall compare with a 
symmetric kernel, using then the radius of the support of J. 

Theorem 4.2. Let J be a compactly supported kernel satisfying (|2.ip . We denote by p the 
size of the support of J : 

p = inf{r > : supp (J) C Br} ■ 
Then the following estimate holds: for any 9 G (0, 1) and T > 0, as R oo, 

sup ^ e-'^«^''^^+°(^''^^) (4.3) 

Notice that J can be asymmetric and have a singularity at the origin. 

Proof. We first use Lemma [2761 to reduce our estimate to the case of symmetric, compactly 
supported kernels. More precisely, using Lemma [2.61 if J* is a symmetric kernel such that 
J ^ J* and supp (J*) = Bp, then for \p\ big enough, ^ L where L* is the Lagrangian 
associated to J*. Taking a look at Theorem 13. H this implies that 

-RIoo{x/R,t/R) = 



— mm t . 

yedBn \ t J 



Now we assume that |a;| < 9R so that in this set \x — y\ ^ {1 — 9)R ^ oo and we shall use 
the behaviour at infinity of L». Lemma 12.51 allows us to wipe out the possible singular part 
near the origin as well as the differential terms of the Hamiltonian. 

Since J* is symmetric, so is L^^^ so that, noting Ll^^{x) = i*(|a;|) we get 

hminf— ^^^iM^i^O. 

yedBR V t 

Then we use the results of Lemma 4.1 and Corollary 4.2] applied to which is symmetric, 
associated to a nonsingular kernel to conclude. □ 



5. Intermediate kernels 



We consider now a general kernel J satisfying (j2.ip . positive everywhere in R^, so that we 
can always write 

J{y) = e-^y^^^y^ , with u;iy) = . 
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We will now make some further assumptions in this section: 

J is C^-smooth for \y\ > 0, 
\DJ{y)\ is bounded on {\y\ > e}, We > 0, 

y ^ \y\u!{y) is supcrlinear and convex, i'^-^) 

3,e(0,l]. hminf 

|yH+oo \y\\Duj{y)\ 

Let us comment these hypotheses: 

(i) The regularity assumption on J (which implies the same regularity for uj) is not crucial 
since by comparison we can deal with less regular kernels, using the results of Section[2l 
(ii) About the convexity of |y|a;(y), it is actually only required for large |y| for the same 
reason: we only care about the tail of J. 

(iii) The case of compactly supported kernels, which would correspond to uj{y) = +00 
outside a ball is treated in Section [H On the other hand, the case of critical kernels 
treated in Section [HI corresponds to Hma;(y) — £ < 00. So, the assumption of super- 
linearity is in between: limw(?/) ~ 00. This is why we speak of intermediate kernels 
here. 

(iv) The superlinearity assumption implies that J automatically satisfies p.4|) since indeed, 
for any /? > 0, we have that uj{y) > (3 for \y\ large enough. 

{v) The convexity and superlinearity assumptions altogether allow us to define the Leg- 



endre transform /C(-) of |y|cij(j/), which will also be super linear and convex, see 2lj : 



/C(p) := sup {p ■ y - \y\uj{y)] . (5.2) 

This function /C(-) will play a big role in estimating the rate of convergence. Notice 
that it is also the Legendre transform of ln(l/J). 
{vi) The "angle" condition on Du says that the gradient cannot take a purely tangential 
position. This is a very weak assumption in this form that allows us to derive a 
minimum behaviour for non-symmetric kernels. 



5.1. Properties of IC 

Thanks to Lemma 12.31 we know that 

L{q)^Pa{q)-DH{po{q)), 

where Poiq) is such that L(q) ~ Pq ■ q — H{po). Thus, a main step consists in finding a lower 
bound for p ■ DH(p). Here is where lC{p) plays an important role: roughly speaking, we will 
see that 

''p-DH{p)= [ p-yef«-l«l'^(^)d?/ ^ e'^^f)" . 

Hence a detailed study of the properties of /C is needed. To this aim, let yo(j)) be the point 
where the sup in (|5.2[) is attained. 

Lemma 5.1. Let us assume that J{y) = c"'^'"'^*') satisfies (|2.ip and (|5.ip . Then the 
function K, is nonnegative and /C(0) = 0. Moreover, for \p\ large enough, the supremum 
in (j5.2p is attained at a unique point yo = yo{p) and \yo{p)\ 00 as \p\ 00. 
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Proof. For p G fixed, let us denote by 

'Ppiy) ■=p-y- l2/k(2/)- 

Then observe that ip{0) = 0, so that IC{p) = sup (pp ^ and of course A^(0) = 0. The 
superlmearity y i— *■ \y\u!{y) imphes that ipp{y) — > —oo as \y\ oo, so that there exists at 
least a maximum point yo{p)- Now, any maximum point satisfies: 

p = D{\yHy)){yo)^^^, 

and by our assumptions, J{y) ^ as |y| oo while \DJ\ remains bounded away from 
y = 0. Thus, as \p\ cxo, necessarily |2/o(p)| +oo. 

We shall now prove that for \p\ large enough, the maximum point yo{p) is unique. Since 
y I— > \y\uj{y) is convex for large \y\, say \y\ > C, then ipp is concave for \y\ > C, independently 
of p. Indeed, this comes from the fact that D'^(pp{y) = ~D^{\y\iLj{y)). Thus we take M > 
large enough so that for any \p\ > Af, any minimum point yoip) satisfies |yo(p)| > C and 
enters the region where ipp is concave. Then, the only case when there may exist several 
maximum points is the case when ipp would be constant on some open set. But since 
Lij{y) — > +00 as \y\ — > oo, this cannot happen for y large. Hence if \p\ is large, the maximum 
is attained at a unique point yo{p). □ 

Lemma 5.2. Under hypothesis (|5.ip . we have 

limmf I II , ^ V- (5.3) 

IpH+oo \p\\yo[P)\ 

Proof. By the definition of yo{p), we have p ~ D{\y\uj{y)){yo) so that 

\p\ ^ uj{yo) + \yo\\Duj{yo)\ 

and 

(yo -p) = |yo|w(yo) + \yo\{yo ■ D^^ivo))- 

This implies, using (|5.1[) that 

(yo-p) ^ v{\yo\^{yo) + \yof\Duj{yo)\) = r]\yo\\p\ , 
hence the result follows taking the liminf. □ 

5.2. Estimating L{q) 

We begin by a technical lemma that will help us in constructing a box where p ■ y — \y\ui{y) 
is close to K,{p) (see below): 

Lemma 5.3. There exists M > such that for any \p\ > M and any ^ G B{yQ{p), l/\p\), 
we have 

\Du;{Om " 4 

with T] defined in (jS.ip . 
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Proof. Notice first that since ^ is at distance at most l/\p\ of (whicli goes at infinity 
as \p\ Qo), tlien all the points ^ G B{yo{p), l/\p\) verify that |^| is big provided \p\ is also 
big. More precisely, using (jS.ip we know there exists a to > such that for any \y\ > m, 

{Duj{y) ■ y) > 377 
\Du;{y)\\y\ 4 ^ " ' 

Now, since |yo(p)I — > cxd as \p\ 00, there exists M > such that if \p\ > M then any 
1^1 S B{yo{p), l/\p\) verifies |^| > \yo\ — > to. Then we may apply the above estimate 
taking y = ^, which gives the result. □ 

Let us now make clear some definitions: 



Definition 5.1. The angle between two vectors a, 6 G is defined as follows: 

ang(a, b) = arccos Y-jni ^ t^^' '^l • 

Moreover, given a vector a G and an angle a G [0,7r/2], we define the positive cone C 
in the direction a with aperture a as follows: 

C+(a, a) = {6 e : ang(a, b) a} . 

Accordingly we define the negative cone as follows: 

C^(a, a) = C^(— a, a) . 



Notice that by considering only apertures ^ a ^ 7r/2 (which is enough for our purpose 
here), we make sure that C~(a, a) D C^(a, a) = {0}. The we have the following lemma: 



Lemma 5.4. Let us consider the cone C*{p) := C (^yo{p),TT/2 — arccos(r//2)) and 
A{p) = {yoip) + z:\z\^^ ze C*{p)} . 

m 

Then for \p\ big enough, we have 

yy e A{p) , p-y- \y\uj{y) ^p-yo - |yok(yo) - 1 • 

Moreover, the volume of A{jp) is bounded by 
for some constant c(iV, 77) > 0. 



Proof. In all the proof we consider at least \p\ > M so that we may apply Lemma [^31 and 
take y G A{p). Then y G B{yo, ^) so that we have [p ■ y) 1^ {p • yo) — 1- 

Moreover, since < rj ^ 1, the aperture of the cone C*{p), arccos(77/2), is strictly positive 
and not equal to tt/2. Hence for \p\ big enough, A{p) C B\y^\ so that \y\ ^ |yo| and it is 
enough to check that uj{y) ^ uj{yo) to get what we want, that is: 

(p-y)- \y\^{y) > {p ■ yo) - |yok(yo) - i ■ 
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To this end, we write 

^(y) - ^(yo) = D(^{0 ■ {y - yo) 

for some ^ G [y,2/o]- The point is that, unless we arc in a symmetric case, we do not have a 
very precise estimate for Duj{^). But we will prove that if y £ C*{p) then necessarily 

DcoiO ■ (y - yo) ^ , (5.4) 
which will be enough to get the desired estimate. 






Figure 1: Cones in yo and ^. 



We shall first show that the image set _Da;(-Bi/|p| (yo(p))) si contained in a cone of aperture 
comparable to arccos(377/4), in the direction yo. More precisely, we claim that for any 
^€B{yo,l/\p\), 

aiig{DLu{^),yo) arccos(37//4) + 0(1 /\p\'^^) , 

see Figure[T] Indeed, we first estimate the angle ang(?/o, for any ^ = yo + z ^ B{yf), l/\p\) 
as follows: 

jyo ■ ^ jyo ■ yo) {yo ■ z) 

\yom |2;o|(|2/o| + 0(l/W)) Iyol(|yol+O(i/H)) 

(recall that |j/o(?')| ^ oo), where the term o(l/|p|) is nonnegative. The expansion, asx 0^, 
arccos(l + x) = 0{x^^^) implies 

ang(yo,C) = o(l/b|i/2) 

and finally, we use the fact that 

ang{Duj{0,y) < ang{DLo{0,0 + ang(e,?/o) arccos(377/4) + o(l/|p|i/2) . 

Thus, for IpI large enough (recall that the arccos function is strictly decreasing on (0, 1)), 

ang(Dw(0,yo) < arccos(77/2) . 

Figure[2]shows the vectorial cone of aperture arccos(377/4) in the direction ^ which is included 
in the vectorial cone of aperture arccos(7y/2) in the direction yo- Hence, if we choose a point 
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y = yo + z such that z G C (yo, 7r/2 — arccos(77/2)) we make sure that ang(_Dw(^), z) ^ 7r/2, 
which yields 

To end the lemma, we only need to mention that A{p) is given by the intersection of the 
ball B{yQ, l/\p\) with a cone placed at yo, of aperture tt/2 — arccos(7]/2). Hence its volume 
is indeed given c{N,i])/\p\^ for some constant c{N,ri) > 0. □ 



Remark 5.5. As rj becomes close to zero, the aperture of the cone C*{p) becomes very 
small: tt/2 — arccos(77/2) ~ 0. So in the limit case 77 = we would not be able to construct a 
non-empty set A{p), at least with this method. Now, the size of A{p) gets small as |p| — > 00 
since it behaves like \p\~^ , but nevertheless this will be sufficient to get a suitable estimate, 
since inside this (small) region, we will get an exponential behaviour. 



Lemma 5.6. Let J{y) = e satisfying (|2.ip and (|5.1|) . Then following estimate 

holds: 

. HpJDI£^%p)) 
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Proof. We start from 

J{\y\>po/2} 



{...)dy+ (...)dy, 

{p-y^O}n{|y|>po/2} J {p-y>0}n{\y\>po/2} 

and the integral over {p ■ y < 0} fl {|y| > /9o/2} is bounded by C\p\ for some C > 0, 
so that it will be negligible (recall that p ■ DH'^^^{p) behaves at least exponentially - see 
Lemma r2.2p . Thus we shall give an estimate of p ■ DH'^^^{p) only in terms of the integral 
over {p ■ y > 0} n {\y\ > po/2}: for \p\ large enough, we have 



p-DH^^^{p)^\ I p-yeP- 

^ J{p-y>0}n{\y\>po/2} 



p-y-lyli^iy) 



dy 



Now let us notice that, for \p\ big enough, A{p) C {p ■ y > 0} Ci {\y\ > po/2}. So we can 
write, using Lemma 15.41 and (j5.ip : 



p-DH'^^^p) ^ I f 

^ Ja(v) 



p.yeP-y-\yMy) 

A{p} 



A{p) 

> ^^HI-4(p)|(|yo(p)|-p)e'=(^)-i 



Hence, taking logarithms and dividing by IC{p): 

ln{p.DH^^^{p)) Hr^\p\/2) , ln|^(p)| , ln(|yo(p)| - l/|p|) , ^ 1 



/C(p) ^ IC{p) 1C{p) 1C{p) ic{p)- 

Recall that lC{p) is supcrlincar. |yo(p)| ~> oo and |-4(p)| ^ c\p\^^ to conclude that 

. ^ \n{p-DH--^{p)) 

□ 

We assume now that J*(y) = c is symmetric in order to get a more precise estimate. 

We denote by L^. and /Ch. the associated Lagrangian and Legendre transform of \y\u)^{y). 
Since in this case /C* is symmetric (and still supcrlincar), wc know that for \p\ large enough, 
/C;7^ : R+ is defined and for any p G JC~^{\z\)^ we have |p| =constant. 

Proposition 5.7. Let us assume that J* is a symmetric kernel satisfying (|2.ip and (|5.ip . 

lim inf ^^4^^ > 1 . (5.6) 

kH- |g||/C.-i(lnM)| 



Proof. In this proof, we drop the ^-subscript for simplicity. We use estimate (|5.5p and, 
since J is symmetric, H°^^{-) is also symmetric. Hence po ■ DH'^^^{po) = \po\\DH'^^^{pq)\ so 
that, using that /C(-) is supcrlincar: 

. ln|po|+ln|Z?i/°^^(po)| .. . ln|Z?i?°-(po)L , 
lim mf — — r = hm mf — - — ^ 1 . 

\q\^^ K-{P0) \qHoo K,[po) 
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In other words, this means that for any e G (0, 1), if \q\ is large enough. 

\n\q\ ^ {1 - e)IC{po{q)) . 

Using again that /C is supcrUncar and symmetric, we have that for any e e (0,1), it is 
enough to take |g| big to get 

bo(g)K|/C-i((l-£)-MnM)|. 
From this we get that for any e G (0, 1): 

"|g||/C-i((l-e)-^lnM)r ' 
and we pass to the limit as e ^ to get the result for L'^^^lq) ~ pq ■ DH°^^{pq): 

liminf , — — — r-rrr 1? 1- 

|,Hoo |g||/C-i(ln|<z|)| ^ 

Finally, we invoke Lemma [2T5l to conclude that the result holds for □ 
5.3. Conclusion 

We are now ready to prove one of the main result of this paper: 

Theorem 5.8. Let J[y) = e~\y\^'^y) be a kernel satisfying ()2.ip and (|5.1|) . Let us consider 
a symmetric kernel J*{y) — e"!^!'^*^^-' such that J J,, and denote by /C* the associated 
Legendre transform of Then the following estimate holds as R —> oo: for any 

61 G (0,1), T> 0, 



sup \u — uii\{x,t) ^ e 



-il-e)RIC-\ln{{l-e)R/t)){l+o{l)) 



Proof. The proof essentially follows the one in the case of compactly supported kernels: 
we first use Lemma 12.61 to reduce our estimate to the case of a symmetric kernel J =^ J* : 

-RLocix/R,t/R) = -R min {t/R)L^^^^^^ 



yedBi \ t/R 

^ -i? min {t/R)LJ^l^^^) 
yeasi ' ' V t/R J 

dist(a;; dBn) 



^ -tL^ 



t 



where i*(|a;|) = i*(x) is the symmetric Lagrangian associated to J». 

Now we assume that |.t| < OR so that in this set dist(x; dBn) ^ (1 — 6)R oo and we use 
the behaviour at infinity of given by (|5.6p to get the result. 

□ 

Corollary 5.9. In particular, if x remains in a bounded set Bm, we can take any 9 > 
and if t remains also bounded we obtain a simpler estimate: 

sup \u-UR\ix,t)^ g-i?K-(lni?)(l+o(l))^ 

|a;[<A/,0<t<l 
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Several remarks are to be made: 

(i) The authors gave in 0] some exphcit estimates, which consists here more or less in 
expliciting the function /C(-). We refer to Section 17.11 for a list of known explicit 
behaviours. 

(a) Even if we are able to prove a lower estimate for asymmetric kernels - Lemma (j5.6p - 
using IC{p) which is not symmetric in general, we are facing a problem: knowing the 
behaviour of p ■ DH{p) is not enough to know the behaviour of each of the vectors, 
unless we make sure that they point more or less in the same direction. And this is not 
clear unless the kernel is "almost" symmetric because of the exponential behaviour. 
This is why we compare with the smallest symmetric kernel above J in order to have 
a more explicit behaviour. 

[in) Even if we were able to derive a bound taking into account the asymmetry, then we 
would have to study the min in Theorem 13. 1[ which is again not obvious unless we 
have an almost symmetric lagrangian. 

(iv) However, see Section [72] for the 1-D case where we can deal with asymmetric kernels, 
since the regions {x > 0} and {a; < 0} are clearly separated. 

6. Critical kernels 

We assume now that J is symmetric and that 

J{y) = e-^y^'^^y^ , where lim uj{y) = £. 

\y\^oo 

Hence J satisfies (|2.5|) with /3o = We want to show that the estimate remains valid 
even if H is not finite everywhere. To this aim, we have to study the Hamilton- Jacobi 
equation more carefully. In this case, the domain of definition of the Hamiltonian H is 
exactly: dom(iJ) = Bi, that is, H = +oo outside Bi. For simplicity, we will first assume 
that £ = Po = 1, the adaptations for other values being straightforward. Then we shall give 
the general result in Theorem 16.71 

6.1. Hamilton- Jacobi equation with nonfinite hamiltonian 

We study the equation ut + H{Du) ~ 0, posed in the cylinder Q := Bi x (0, oo) (although 
most of the results of this section would hold also for more general cylinders). Here, we 
assume that the hamiltonian H is infinite in the complement of which is the main 
difficulty. Following we begin by constructing a new equation which is equivalent in 
the viscosity sense to ut + H{Du) = 0, the main interest being that it allows us to prove 
comparison and analyze the initial trace of solutions. 

On the parabolic boundary 

dpQ = (dBi X [0, oo)) U (Bi x {t = 0}) , 

we impose a continuous boundary condition / in the viscous sense. More precisely, we 
consider the following problem: 

'ut+H{Du) = inBix(0,oo), 

< u{x,t) = f{x,t) on 9Si X [0, oo) , (6.1) 
u{x,0) = f{x,0) inSi. 
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Definition 6.1. Given / G C{dpQ), wc say that an upper semi-continuous function u is 
a viscosity subsolution of (|6.ip if for any smooth function ip such that u — (p reaches a 
maximum at (xo,to) we have: 

{xo,to)eQ =^ tpt + H{Dtp) 1^0 at {xo,ta) 
{xq, to) e dpQ min (^tpt + H{Dtp) ; m - ^ at (xq, to) . 

The same definition holds with reversed inequahties (and min instead of max) for an upper 
semi-continuous viscosity subsolution. And finally: 

Definition 6.2. A locally bounded function m : Q ^ M is a viscosity solution of (|6.ip if its 
upper semi-continuous envelope is a supersolution and its lower semi-continuous envelope is 
a subsolution of (|6.ip . 



Let us mention that in the case when H is finite everywhere, solutions take on the initial 
data in a classical way. But since here some data may not be compatible with the fact that 
dom(iJ) = Bi, this implies that a boundary/initial layer appears, and this is precisely the 
phenomenon we are facing. In order to understand this layer, we need first to reinterpret 
the equation with a new Hamiltonian: 



Proposition 6.1. Subsolutions and supersolutions of (|6.ip are also suhsolutions and su- 
persolutions (in the viscous sense) of the equation: 



max l^ut + H{Du) ; \Du\ - 1 j = . (6.2) 
with the same data on the parabolic boundary. 

Proof. Let u be a viscosity subsolution and consider a smooth test function such that 
u — if has a maximum at (xp, to). We assume for simplicity that (.tq, to) G Q (the argument 
being similar if it is a boundary point) . Then since by definition 

dMxo,to) + H{Dip){xo,to) ^ , 

we have necessarily H{Dip) < +oo, so that \Dip\ ^ 1 and thus u satisfies (in the viscous 
sense) also the inequation 



max [dtu + H{Du) ; \Du\ - 1 j < . 

Now if u is a supersolution and ip is such that v ~ ip has a minimum at (xq, to): then 

dtip{xo,to) + H{Dip){xo,to) ^ 

implies that 

v[iax{dtip{xoM) + H{Dip){xo,to)\ |D(p(xo, to)| - l) ^0, 
which implies that w is a supersolution of (|6.2p . □ 

Proposition 6.2. Let u be a viscosity subsolution and v be a viscosity supersolution of 

max [ut + H{Du) ; \Du\ - 1^ = . 
with u ^ V on the parabolic boundary dpQ. Then m ^ w. 
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Proof. Formally speaking, if one fixes fi G (0,1) and considers a maximum point (x'o,io) 
of fiu ~ V, then if iZJuK^o, to) ^ 1, since by definition fiDu ~ Dv at (xq, to), it comes that 
\Du\(xo-,tQ) 1/ fi > 1 which is impossible since m is a viscosity subsolution. So we have 
both \Dv\, \Du\ < 1 at (xq, to) and we do the comparison as always, using standard viscosity 
techniques. Now let us be more precise. 

We fix /i e (0, 1); and T > and consider a point (xq, yo, to, sq) G -Bi x [0, T]^ where 

r ^ , ^ \x-y\^ \t-s\^ C 
<P : (x, y, t, s) ^ nu(x, t) - v{y, s) — 



{jf^ + + Hi^,-'p) ; VI - l) ^ . (6.3) 



/32 T-t 

reaches its maximum. We assume that it is an interior point otherwise, using the boundary 
values one obtains immediately <i> ^ in i?^ x (0,T)^ which is what we want. 

Fixing one variable, since {x,t) ^{x,yo,t, sq) reaches its maximum at (xo,to) one may 
consider the following test function for u at {xq, to): 

{X, t) - (.(,„, .o) + + + J— J ■ 

Indeed, if we denote by p :— 2{xo — yo)/£^, q '■= 2(to ~ so)//9^i it comes 

,(r^2 

On the other hand, for v we use at (yo, so) the test function 

, . ,N . , N No - y? \to - 

92(3, t) := fiu{xo,to) ^2 ^2 — ' 

which leads to 

max (^q + H{p) ; |p| - ^ . (6.4) 

If we assume that \p\ ^ 1 then ^ > 1 which is impossible from (|6.3p . So, both \p\ 

and are less than 1 and then the proof follows standard arguments of viscosity solu- 

tions: we can combine (|6.3p and (j6.4p . getting rid of the max which gives (after multiplying 
the first inequality by fi): 

^ fiH{p-^p)-H{p)^0. (6.5) 



(T-t)2 



We claim that /i(^) := iiH{iJL^'^p) — H{p) ^ for any p G and ^ G (0, 1), which leads to 
a contradiction with (j6.5p . so that an interior maximum of $ is impossible. Hence, $ ^ 
in (Bi)^ X (0,T)^ and since P,e,C,T, fi are arbitrary, we finally conclude that u ^ v in 
Bi X (0,oo). 

To end the proof, let us check the claim: using the convexity of H, one gets 

h'{^i) = ((m" V) • DH{p-'p) - H{ii-^p)) ^ , 
and since h{l) — 0, we see that /i ^ for /i G (0, 1). □ 

Proposition 6.3. Let u he a viscosity solution of max (^Ut + H{Du) ; \Du\ — 1^ =0, and 
u = f on dpQ. Then the initial trace of u, tt(0) verifies: 

max(u(0)-/; |i:'u(0)| -l) =0. (6.6) 
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Proof. Let u be a subsolution of the equation with boundary data /. Then at i = 0, 
min max('Ut + H{Du), \Du\ - 1) ; u - s$ . 

We take a test funetion (p{x, t) ~ CA +\x ~ xop/e^ such that 

max((pt + H{Dip), \Dlp\ - 1) ^ . 
This is always possible if e is small enough, so that this implies ii(0) ^ /. Thus: 



(u{Q)~f - \Dum 



I K 0. 



Now, we consider a supcrsolution u and take Lp{x) such that f — (p has a minimum at xq. 
Then for any C > 0, 

ip : {x, t) (-^ u{x, 0) - ip{x) + Ct 
has a maximum aX t = Q, x — xq- Using tjj as test-function, we obtain 

mayi(^mayi{-C + H{Dip{0)),\Dip{0) \ - 1) ; u - ^0. 

For C big enough, we have — C + H{Dip{0)) < so that there remain two possibilities: 
Case 1: \Dlp{0) - 1| < 0, which implies m - / ^ so that: 

max((^(0)-/; |Z?(^(0)| - l) ^0. 

Case 2: \Dlp{0) 0, which implies 

max((^(0)-/; \D^m - l) >0. 

So if u is a solution, both inequalities give the equality. □ 

Remark 6.4. Equation (|6.6p can be understood as an obstacle problem: both |£'u(0)| ^ 1, 
u{0) ^ f and (m(0) - f){\Dum - 1) = 0. 

6.2. Back to the estimate of I 

As in Section m we define for any ^ > 0: 

I^:=--ln{e-^^ + l), 7^ := lim sup , := lim inf . (6.7) 

Taking liminf and limsup and a test function, we see that for any A > 0, and are 
respectively sub and supersolutions of 

max (ut + H (Du) ; \Du\-l^ =0, I^{0)^A. 

. 

Consequently, using (|6.2p . we obtain that I = 1 , so that as i? — > -t-oo, all the sequence 
/j| converges to the unique solution /"^ of the problem. Then as A +oo, ^ I which 
satisfies the equation with /(0,x) = dist(x; 

Now we have to identify the limit /, and to do so we have to study some properties of this 
specific Lagrangian. 
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Lemma 6.5. The Lagrangian L satisfies the following properties: 
\DL{q)\ < 1 and L{q) ~ \q\ as \q\ ^ oo . 

Proof. Since by definition, 

L{q) = sup {p-q- H{p)} = snp {p ■ q - H{p)} , 

pGK™ |p|<l 

and H{p) +00 as \p\ 1, this implies that the sup is attained at some |po('?)| < 1- 
On the other hand, a simple calculus shows that DL{q) = po{q), so that indeed, for any q, 
\DL{q)\ < 1. This also implies a first basic estimate: L{q) ^ To get the equivalent, we 
first bound H{p). Let us first notice that 

H{p) = H-'%p) + fip) = if-^(p) + 0(1) , 

indeed the singular and differential parts of the hamiltonian remain bounded in the set 
{|p| < 1}, as well as Df{p). Now, for any \p\ < 1, 

H{p) ^ f{p)+ f e (l^-l-i)!^! dy ^ f{p)+c(N) H e (M-^K^-^ dr = / (p) + ^i^^^K^^^-iH . 

Jr^ Jo U ~ \P\) 

and this sup is attained for pq = po{q) satisfying the equation: 

q = \p\c(N){N - 1)!(1 - bl)""^"'^ + Dfip) . 
Thus, as \q\ — > 00, necessarily |po| — > 1 and since po and q point in the same direction, 

L{q) > sup {p • (? - '^''ff''^^ - fip)} ~ Po{q) • g ~ kl • 
|p|<i U ^ \P\) 

Since we have seen that L(q) ^ \q\, we conclude that 

L{q)^\q\ as |g| -> 00 . 

□ 

Since J is assumed to be symmetric, so are H and L and we write L{\x\) ~ L{x). 

Proposition 6.6. The solution of ()6.1|) with initial data 1(0, x) ~ dist{x;dBi) and 1 = 
on the boundary is given by the Lax- Oleinik formula: 

dist(a;; dBi)^ 



/(.,t) = tL(=:^). 



Proof. Since \DL\ < 1, then \DI\ < 1 so that the compatibility condition is always fulfilled 
and the equation holds everywhere. Now we take a look at the initial data. Since L{q) ~ \q\, 
this implies: 

tL[^^^^^)^dist{x;dB,)ast^O, 
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so that indeed the Lax-Oleinik formula gives a solution. Since the viscosity solution is 
unique, this ends the proof. □ 

The reader will easily check that ii £ ~ /3o ^ 1, then all the results of this section remain 
valid and then ^ /3o|9| a-s \q\ — > oo. Moreover. Lemma 12.51 and Lemma 12.61 are also 
valid in the present case since |p| remains bounded: as we have seen, = H'^^^ijj) + 0(1) 

as \p\ = Po- 

Hence we may write down a more general result for possibly non-symmetric and singular 
kernels: 

Theorem 6.7. Let J be a kernel satisfying p.ip and p.5[) . In particular, J can be asym- 
metric and have a singularity at the origin. Then for any 9 E (0, 1), we have the following 
estimate as R oo: for any 9 € (0,1) T > 0, 

sup |u-Mi?J(a;,i) ^ e-(i-^)*«+°(^) (6.8) 

Proof. We skip the details since this is the same as for Proposition 14.21 we first reduce 
the estimate to symmetric kernel by comparison, putting a symmetric kernel above J with 
the same (3q in (|2.5p and then we wipe out the possible singularities and differential terms. 
Actually, the proof is even simpler since since those terms remain bounded in the set {\p\ < 
Po}, hence only the exponential part of the Hamiltonian plays a role in the estimate. □ 

Remark 6.8. As Pq — > 0, the estimate gets worse. Indeed, this means that the kernel tends 
to behave slower than an exponential and we are facing a problem of fat tails (like a power 
decay), that this method cannot handle. 

7. Further results and comments 
7.1. Explicit bounds 

In some particular cases, we already gave concrete estimates in @, by studying directly the 
Hamiltonian H. Computing the function IC~^, we recover here these estimates under the 
following form: 

sup \u-un\{x,t) s; e-«'c"('"«)(i+°(i)) 

|2;|<A/,0<t<l 

Table \7A\ collects some known asymptotic behaviors of IC~^ and sup(u — uji). Notice that 
since K, is superlinear (except in the critical case), then IC~^{z) is defined for z > big 
enough. Most of the calculations are straightforward, we only sketch the case J{y) = e ° " : 
in this case, IC{y) \y\ In |?;| and if we define z = \y\ In \y\ we get Inz = In \y\ + In In |y| ~ 
ln|j/|. Hence z |y|lnz. which implies \y\ ~ z/lnz. Let us also mention that in the case 
J{y) = e ""Ifl, IC~^{z) = a in the sense of graphs. 

Remark 7.1. As we have seen, the presence of a singularity at the origin does not modify 
the behaviour of IC~^, so by instance, if we multiply by \y\~'^ any of the previous kernels we 
obtain the same estimate. 
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Kernel 



sup \u — ub\ 



Jiy) = ^{\yKp}iy) 

J(y)=e-'"' 
J{y) = c-^y^", a > 1 

J{y) = c-lfl'"l?^l 
J{y) = e-°'\v\^ a > 



zp 

In z 
a 



g-flln(i?,)/p 
g -R{R\nR} 
-i?,(lni?,)(°-i)''° 
g -_Rlnln_R 
„ -aR 



Table 1: Behaviour of /C ^ 



7.2. Non-symmetric kernels in one space dimension 

As we have seen, it is in general difficult to give an explicit behaviour of sup \ u — uii\ for 
non-symetric kernels, unless we compare with a symetric one. But in the case A'^ = 1, the 
regions where {x > 0} and {x < 0} are clearly seperated so that more precise estimates 
can be given according to the tails of J at — cx) and +00 which can be different. We shall 
just illustrate this in an explicit case which concerns the most extreme cases we can cover: 
on the one side we have a compactly supported kernel, while on the other side, we have an 
exponential decay. 

Proposition 7.2. Let us assume that N ^ 1 and J is given by: 

Jiy) — ^e-l2^ll|y<o} + ^l{a^o}- 
Then the associated Lagrangian satisfies: 

L{q) ^ '/lug', and L{q) ~ q. 

q — ^ + 00 q — ^ — 00 



Proof. A straightforward calculus shows that the Hamiltonian is defined in the region 
{p > —1} and that: 

eP 1 1 

Hip) = — - — + , - 1 . 

2p 2p 2(39+1) 

For q > 0, we calculate the Lagragian as follows: 

L{q) = sup {pq - H{p)} ^ sup{pq - H{p)}. 

P>-1 p>0 

Indeed, if q > 0, then pq > only for p > and we know the sup is nonncgativc so that 
it has to be attained for p > 0. With this remark, the estimate is just the same as for the 
case when J{y) = ^^{y^o}- The same remark holds when q < 0: the sup is attained in the 
region — 1 < p < and the behaviour is given by then exponential decay of J. □ 

Then we are able to use Theorcm l3.1l in a more precise way (we consider x fixed in a bounded 
domain for simplicity): 
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Proposition 7.3. Let J be defined as above. Then 

sup \u-uii\{x,t) 

0<2;<Af,0<t<l 

sup \u-uii\{x,t) ^ 

-j\f<a;<0,0<t<l 

Proof. We just come back to the expression of loo- 

Jo,(x/i?,</i?)= min 

So, if 1 < X < M, whether y = —R or y = +R, we always have Rx ~ y ^ {M — l)R +00. 
Hence the min of L is attained for y — ~R and we recover the behaviour of L{q) for 
q = R{M + 1) — > 00, that is, a i?lni? behaviour. On the other hand, if — M < < — 1, the 
min is attained for y = —R and we get the Unear behaviour of for gr ^ — od. □ 



g -filn_R(l+o(l)) 









r \ 







75 U U U , , , U U u 

-20 -16 -10 -5 5 10 15 20 



Figure 3: Convergence for non-symmetric kernel. 



In Figure [3] we plot for uo = 1 the approximations ur for R ~ 10, 15, 20 and non symmetric 
J as in Proposition 17.21 This illustrates the different rate of convergence whether a; < or 
X > 0. 



7.3. KPP-type results 

In this section, we briefly explain how our results allow to treat a non-local version of 
the KPP-problem (Kolmogorov-Petrovskii-Piskounov) associated to equation (|l.ip with the 
classical monostable u{\ — M)-term. For simplicity we shall just explain this on the following 
equation: 

dtu- [J *u-u) =u{l-u) in K^x[G,T], 

with a continuous initial data uq{x) = g{x), 5j g ^ 1. Existence of solutions with initial 
data ^ .9 ^ 1 may be obtained for instance by Perron's method. 

The interested reader will find further references about this equation and traveling waves in 
the works of J. Coville and L. Dupaigne [l6l |. 
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Now, in order to study convergence of u to the equilibrium states and 1 for large x and t, 
the following scaling is widely used: 



It turns out here that formally, satisfies the equation 



dtu' - i ( Je * - u") = ^ — ^ in X [0, T] , 

e £ 

with Je{x) ~ e'^J{x/e), so that one may use exactly the same method as was used in 
Section [3] with e playing the role of l/i? ^ 0. 

We may thus combine the techniques of [1] with the ones we used in Section [3] to handle 
the convergence of the non-local term, which give us some estimates at which convergence 
to the states u = and u = I occur. In fact, once we know how to deal with the non-local 
terms of the equation, the rest of the proof only follows this is why we only sketch a 
proof below. Here we denote by H the following Hamiltonian: 

H{p) := / (eP y -l)j{y)dy, 
where it is assumed that J G C(R^) and satisfies (j2.5p for some Po > 0. 

Theorem 7.4. Let Go = {9{x) = 0} and Gi = {g{x) = 1}. Then the following results 
hold: Let Lq ~ — eln(u^). Then Iq Iq, the solution of the variational inequality: 

'di. 



u-i(^-^+H{DIo) + l, Iq^ =0 m K'^x(0,oo), 



and 



loix,0) 



+00 if X G Go , 

ifxem^\Go. 



(a) Let Lf := — eln(l — u'^). Then If Ii locally uniformly in the set C = {ln{x,t) = O}, 
where 

— i + H(Dh) + 1 = m int C 
ot 

/i(a;,t) = on dC H {t > 0} 

Ji{x,0) = +cx)l{Gi}(x). 

Sketch of proof. For (i), it turns out that u"^ satisfies the equation 

dtu' ~-{J,*u' - u') = """^^'"'^ in X [0, T] , 

e £ 

with initial data u^{x^{)) ~ g{x/£), where the rescaled kernel is Je(x) = £^J(x/e). 

Then we make the log-transform 

Iq^{x, t) = — eln {if{x, t) + cxp ( — — ) j 

which satisfies 

dJl^^+ I (e-HA'-(^+-Wo-(^)}-l)j(y)dy = ^^^^^TT 

" 7 V ) u-+exp(-A/£) 



35 



Passage to the limit in the left-hand side is done exactly as in Section [3] while handling the 
right-hand side follows exactly from KPP classical techniques: first notice that by construc- 
tion, Iq'^ ^ 0. Then, if in the limit I^ix, t) = lim^^o I^'^ ^^^^ means that ^ so 
that clearly as e ^ 0, the right-hand side converges to 1. 

For (ii), we set similarly 

/i'^(a;,<) = -eln(^l - u''(a;, t) + exp ( - 

The first step consists in proving that limsup^^Q /^''^ and if := liminf£^o ^i'^ 

respectively sub- and super-solutions of the variational inequality: 

max (dtl + H{DI) + 1,1 - ^/'^) in x (0, oo) , 

/(x,0) = A-t{G,}{x), 

where = 1 if {x,t) G C and otherwise. Direct comparison between if and if 

cannot be derived here since no information on the regularity of C is available. The final 
result then follows from a representation formula for if. We refer to 0] for the details. □ 

Then, coming back to the original variables, one can obtain explicit exponential convergence 
rates, which follows from our study of the asymptotic behaviour of the Lagrangian associated 
to H (as was done in Sections SI [51 . 

7.^. Relation with optimal existence results 

We would like to add another final comment on a related subject. As was said, by estimating 
sup^^ \u — ub\{x,T), we are measuring the total amount of processes that can escape the 
box Br between t = and t = T . Another way of understanding this is that we are somehow 
estimating the Green kernel associated to the equation. 

Thus in [lo| . the authors together with R. Ferreira are deriving similar estimates but in 
the context of optimal initial data, which is also a way of measuring the behaviour of the 
kernel at infinity. Hence it is not so surprising that similar estimates appear, even if they 
are obtained through a totally different method. 

For instance, it turns out that if J is compactly supported, the optimal class of existence 
for = J * u — u in consists of initial data satisfying: 

\uo{x)\ ^ Ccl-I"^l-I, 

hence we recover a -R In i? estimate for the Green function, typical of compactly supported 
kernels. We refer to for more results in this direction. 
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